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ABSTRACT 

In this work multi phase flow and erosion analysis 
were done via simulations & CFD application for 
standard diesel and two alternative bio fuels, FAME 
and DME, inside different nozzle models and with 
various boundary conditions. Nozzle model consists 
of narrow channel with sharp (type I) or rounded 
(type Y) inlet section, with or without downstream 
placed target, so there was a total of four different 
model geometries. Simulation results showed that 
cavitations was present in almost all cases and that 
clear difference between three observed fuels can be 
seen. Mass flow in channel type I was lower than one 
in channel type Y. When comparing three observed 
fuels, it was noticed that DME fuel usually had 
highest velocity, but lowest mass flow rate. Contrary 
to DME, FAME fuel showed highest mass flow rate 
despite lowest velocity. When comparing simulation 
results and physical properties of observed fuels, it 
can be concluded that density is a leading term in 
determining mass flow rate. Also, erosion model 
predicts more intensive MDPR value near narrow 
channel exit. 
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INTRODUCTION 

Alternative fuels are lately becoming more and more 
interesting due to a fact that they don't contribute 
increasing of carbon in atmosphere in Earth's 
atmosphere carbon-cycle. In Europe and U.S. fuels 
gained from rapeseed (RME) and soybean (SOME), 
together called FAME fuels or biodiesel, are used as 
an alternative fuel or as a compound in standard 
diesel/biodiesel mixture. Another biofuel is very 
interesting, DME, which can be gained from all sorts 


of different sources. In Europe there are currently 
several active norms that are regulating amount of 
biodiesel in mineral diesel (EN 590 allow up to 5%), 
as well as properties of FAME fuels (EN 14241). 
Diesel fuel injection equipment manufacturers 
brought common statement in which they support the 
development of compression ignition alternative fuels. 

Diversity of physical properties between mentioned 
fuels causes their different flow characteristics inside 
fuel nozzles. In this work flow simulations with 
analyses were performed in different nozzle models, 
for every inlet/outlet pressure drop and three different 
fuels. Calculations were done via CFD. Analysis is 
performed by following criteria: volume fraction 
distribution due to cavitations, achieved mass flow 
rate and velocity profile in narrow channel section. 
Additional figures are presented: pressure, absolute 
velocity and turbulence ki netic energy distribution. 
Also, erosion modelling was done. 

For simulations is used one of the commercially 
available CFD application, AVL's CFD Workflow 
Manager with FIRE solver (hereafter 
CFDWM/FIRE). In CFDWM/FIRE one can perform 
simulations of multiphase flows via Multiphase 
module with Erosion modelling included. 

LITERATURE REVIEW: 

There are some papers which have been studied and 
referred on my work. 

Junmei Shi et al. CFD investigation of fuel 
properties effect on cavitating flow in generic nozzle 
geometries Fuel saturation vapour pressure variation 
up to 0.8 bar has know essential effect on cavitation 
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this is due to the very high pressure gradient caused 
by the flow acceleration in the nozzle. 

Takenaka et al. CFD analysis of nozzle cavitations 
Investigation of cavitation phenomena showed that for 
nearly all the duration of the injection process the 
nozzle injector hole is surrounded by cavitation films. 

Berchiche et al For interfacial mass exchange Linear 
Cavitation Model was used and for interfacial 
momentum exchange Cavitation Drag Model was 
used. Both interfacial exchange models imply two 
additional transport equations: Bubble Number 
Density and Interfacial Area equation. These 
equations bring up additional closure coefficients of 
mathematical model. 

Semelsberger et al FAME fuel achieves highest mass 
flow rates in all cases. DME has lowest mass flow 
rates. Generally, mass flow rates in Channel Y cases 
are higher than ones in Channel I cases for same 
boundary conditions. Both interfacial exchange 
models imply two additional transport equations: 
Bubble Number Density and Interfacial Area 
equation. These equations bring up additional closure 
coefficients of mathematical model. 

Ornella Chiavola et al On a modified VCO nozzle 
layout for diesel common rail injectors under actual 
needle displacement FAME fuel achieves highest 
mass flow rates in all cases. DME has lowest mass 
flow rates. Generally, mass flow rates in Channel Y 
cases are higher than ones in Channel I cases for same 
boundary conditions. 

METHOD USED 

> In this work quasi-stationary, inner, non- 
compressible, viscous, turbulent and two- phase 
(generally multiphase) type of flow is assumed, so 
simulation is set-up considering these 
assumptions. Two-phase flow implies bubble 
number density and interfacial area between liquid 
and gas phase, so additional transport equations 
for modelling cavitation should be active in order 
to describe phase change. 

> All conservation equations can be written in 
generic form Turbulence model is a part of a Eddy 
Viscosity Models (EVM) with two equations 
EVM assumes direct analogy between molecular 
and turbulent momentum. Based on Boussinesq 
assumption. 

> Erosion is in Multiphase module modelled within 


two model quantities Erosion Incubation Time 
Mean Depth Penetration Rate (MDPR) 
Dimensions that are strictly defined: nozzle model 
thickness, narrow channel, target and height (y 
coordinate). Length (x coordinate) of inlet and 
outlet zone is not strictly defined and should be 
selected to give developed flow on inlet and outlet 
zone Calculation stability is also affected with 
these parameter. 

> Because of overall of four different nozzle model 
types, and computational domains respectively, 
hereafter computational domains will be referred 
as: 

• Channel I —> channel type I without target 

• Channel Y —> channel type Y without target 

• Target I —> channel type I with downstream 
placed target 

• Target Y —■» channel type Y with downstream 
placed target 

SELECTIONS AND BOUNDARY CONDITIONS 

All domain meshes have same number and type of 
selections. Boundary conditions must be defined for 
every phase separately, except pressure which is 
coupled variable and is defined only for continuous 
phase. Also, total of phase volume fractions must be 
equal to one. Since driving force for flow quantities is 
pressure difference between inlet and outlet selection, 
in order to avoid too long mesh on inlet zone, flow 
direction on inlet selection is determined. Boundary 
turbulence parameters were pre-defined by AVL. 

CONCLUSIONS 

CFDWM/FIRE application can give us valuable data 
regarding multiphase cavitating flow inside nozzles, it 
is shown that mathematical model is valid for its 
purpose and that nozzle flow for different fuels can be 
done, when comparing different fuels simulation 
results, some obvious differences can be observed. 
Most significant conclusion regarding to comparison 
between fuels is that in mass flow rate equation, fluid 
density is more significant term than molecular 
viscosity. This is concluded based on fact that in 
almost all cases DME fuel had highest velocity but 
lowest mass flow rate. If we combine that with fluid 
properties above conclusion can be derived. 
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